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Abstract 

The non-diagonal correlators of vector and scalar currents are considered at three- 
loop order in QCD. The full mass dependence is computed in the case where one of 
^ • the quarks is massless and the other one carries mass M. We exploit the decoupling 

^ relations between the full theory and the heavy quark effective theory (HQET) in 

order to obtain the logarithmic parts of the leading threshold terms. With the help 
of conformal mapping and Fade approximation numerical estimates for the non- 
logarithmic terms are extracted which in turn lead to a prediction of the correlator 
in HQET at order . As applications of the vector and scalar correlator we consider 
the single-top-quark production via the process qq tb and the decay rate of a 
charged Higgs boson into hadrons, respectively. In both cases the computed NLO 
corrections are shown to be numerically much less important than the leading ones. 
On the contrary, the NLO order QCD corrections to the HQET sum rule for the 
leptonic decay rate of a heavy-light meson proves to be comparable to the leading 
one. 
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1 Introduction 



In QCD the correlator of two currents is very often a central object from which important 
physical consequences can be deduced. In particular, in case the coupling of the currents 
is diagonal physical observables like e^e~ annihilation into hadrons and the decay of the 
Z boson are obtained from the vector and axial-vector current correlators. Furthermore, 
total decay rates of CP even or CP odd Higgs bosons can be obtained by considering the 
scalar and pseudo-scalar current densities, respectively. For these cases the results are 
available up to order a1 taking into account the full quark mass dependence both for the 
non-singlet |jl|] and singlet correlators. 

In this work we want to extend the techniques developed in 0] to the situation where 
the coupling of the external currents to quarks is non-diagonal. In particular we compute 
the three-loop correlators of currents which couple to two different quark flavours. To 
be precise we want to consider the case where one of the quarks is massless and the 
other carries mass M. In this limit the vector (scalar) and axial-vector (pseudo-scalar) 
correlators coincide. A short account of our results has been presented in 0. In this 
paper the details of the calculation are given and explicit expressions for intermediate 
results are provided which might be useful for other applications. 

The main aim of this work is to obtain results which are valid for arbitray values of 
the quark mass, M, respectively the ratio g^/M^ where q is the external momentum of 
the correlator. Very often the knowledge of the full mass dependence is not necessary. 
E.g., in high energy experiments where the center-of-mass energy, a/s, is much larger 
than the mass of the quarks, the latter can often be neglected or an expansion in M/ y/s 
is sufficient to describe the data. However, there are situations where the full dependence 
on M and y/s is required. One can, e.g., think on the high precision which meanwhile has 
been reached at LEP (CERN), SLC (SLAC) or TEVATRON (Fermilab) or on situations 
where the center-of-mass energy is of the same order of magnitude as the quark masses. 
In particular for threshold phenomena the masses are important and the full dependence 
is desirable. 

Concerning the physical applications of the non-diagonal vector and axial- vector cor- 
relators we have in mind properties connected to the W boson. With the results of this 
paper at hand a certain (gauge invariant) class of corrections to the Drell-Yan process 
becomes available. In particular we have in mind the production of a quark pair through 
the decay of a virtual W boson generated in pp collisions. The absorptive part of the con- 
sidered correlator is directly related to the decay width of the (highly virtual) W bosons 
into quark pairs and gluons. Of particular interest in this connection is the single-top- 
quark production via the process qq tb. The imaginary part of the transversal W boson 
polarization function constitutes a gauge invariant and finite contribution of 0{a1). 

As an application of the scalar and pseudo-scalar current correlator we want to mention 
the decay of a charged Higgs boson which occurs in extensions of the Standard Model 
(SM). The corrections provided in this paper describe the total decay rate into a massive 
and a massless quark. 

Another important application of the non-diagonal vector and scalar current correlator 
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is connected to the corresponding meson decay constant. Within HQET it is related to 
the imaginary part of the effective current correlator, i.e. the spectral density, evaluated 
near threshold. As the current in the effective theory does not depend on the 7-matrix 
structure it can be derived both from the vector and the scalar correlator of the full 
theory which are evaluated in this work. The corresponding analysis at order as has 
been performed in Ref. Q with the conclusion that the determination of the B meson 
decay constant suffers from large pertubative corrections. We show in this paper that the 
corrections of order are also sizeable. 

The paper is organized as follows. In Section ^ we provide useful definitions and 
describe the method we use for the computation. The connection between the effective 
and full theory is established in Section |^. In Section ^ we present explicit results for 
the moments of the vector and scalar current correlator and discuss the results for the 
corresponding spectral densities in the full theory. Afterwards, in Section ^, the procedure 
is discussed which allows for the extraction of threshold information from our results. This 
can be translated to the effective theory which provides a result for the spectral density 
in the effective theory up to order a^. Finally, in Section ^, we discuss some physical 
applications and present our conclusions. The Appendix contains the analytical results 
of the low- and high-energy moments for the vector and scalar correlator. 

2 Conformal mapping and Pade approximation 

In this Section we want to describe the method we use for the computation of the current 
correlators in full QCD. Let us start with some definitions. 

In the vector case the polarization function is defined through 



{-q^9,u + q,qu)Ti^{q^) + q,quTil{q^) = ijdxe'^%Q\Tjl{x)j:\m). (1) 



with j'^ = (m(/i)/M)?/'i?/'2 where m(/i) is the MS and M the on-shell mass of tp2- ^1 we 
consider to be massless. Thus we will identify i/ji with q, the massless quark, and ip2 
with Q which is supposed to be a heavy quark of mass M. This will become relevant 
later on when we consider the effective theory. Throughout this paper we consider anti- 
commuting 75 which is justified as for ipi 7^ 1^2 only non-singlet diagrams contribute. As a 
consequence the axial-vector (pseudo-scalar) correlator coincides with the vector (scalar) 
one. 

It is convenient to introduce the dimensionless variable 



with = ipi'~ffj,ip2- Only the transversal part n*'(g^) will be considered in the following. 
The definition of the scalar polarization function reads 




(2) 



z 




(3) 
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where M refers to the pole mass. For the overall normalization of Ii^{q^) {5 = v, s) we 
adopt the QED-like conditions n^(0) = 0. 

The physical observables R{s) are related to n(g2) through 



where the use of the variables 



X 



127r Im 
Stt Im 



M 



s + ie] 



Y[\q' = s + it) 



(4) 
(5) 

(6) 



turns out to be useful to describe the high energy and threshold region, respectively. 
The expansion of 11'' in terms of as reads {5 = v,s) 



n(o).' 



a 



(rif) 



TT 



TT 



(7) 



It is convenient to further decompose the three-loop term according to the colour structure 



n(2). 



r{2),5 



(2),5 



r(2),<5 
''FH ' 



where analogous formulae hold for R{s) and also for the corresponding quantities in 
the effective theory which we will define below. Cp = {N^ — l)/{2Nc) and Ca = 
are the eigenvalues of the quadratic Casimir operator of the fundamental and adjoint 
representation, respectively, and T = 1/2 is the index of the fundamental representation. 

fS) S ('2) S 

In Eq. (^) Upp corresponds to the abelian part already present in QED whereas n^j( 
represents the non-abelian structure. The remaining two structures correspond to the 
fermionic contributions where ni counts the number of massless quarks and Uf = ni + 1 
is the total number of active quark flavours. 

For later convenience we want to list the exact expressions for the Born results and 
the corrections of order a, 



J' L_l' M' L_l' 



10 . In the vector case we have 



2 

3N, 
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(l -x^)^ (2 + x^) 
1 



Sx^ 2x^ Sx'' 1 

— + — + — + -x^{-5 - 4x2 + 5x^) ln(x2) 
2 6 b 6 



- -(1 - x2)2(4 + 5x2) _ ^2^ 

3 



2 
3' 



- ^(l-x2)2(2 + x2) 



- In 



X 



ln(l -x2) + 2Li2 



X 



(9) 
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Correspondingly, the scalar current correlators read 



-^(1 - x^)l - (3 - ^ 2x4) ln(x2) 



2(1 -x' 
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+ (5 - 2x2) ln(l - x^ 



In 



X 



1-X2 



ln(l -x2) + 2Li2 



X 



1 - x^ 



■ (10) 



A complete analytical computation of 11'^ (g2) at three- loop order or its imaginary part 
is currently not feasible. The method we use for the computation of the diagrams has 
successfully been applied to several physical quantities (see, e.g., Refs. |jTT], [l|, 0, |12|)- It 
allows for the computation of a semi-numerical approximation for n(g2). The aim is the 
reconstruction of the function n(g2) from the knowledge of some moments for z — > and 
z — >■ — oo and additional partial information about the behaviour of R{s) for s — > M'^. 
For convenience we briefly summarize the different steps which have to be performed for 
the individual pieces of Eq. (|^) and (||). In the following we generically write n(g2). 

1. Compute as many moments as possible for small and large z. In our case the 
expansion for z ^ reduces to a simple Taylor series of the Feynman diagrams in the 
external momentum. For z ^ —oo, however, the rules of asymptotic expansion 
have to be applied. Thus in this limit besides the power corrections there are 
logarithmic terms in z. 

2. The information known about the imaginary part R{s) for s — > M'^ has to be 
transformed to a function 11*''^ (g2). Afterwards 11*^'" (g2) has to be expanded in the 
limits z —>■ and z ^ —oo and the moments have to be subtracted from the ones 
of n(g2). It is important to construct 11*''^ (g2) in such a way that its expansion 
for 2; — i> is analytical. In this way the information about the logarithmic part of 
n*'"'(g2) can be incorporated. 

3. Construct a function n'°^(g2) in such a way that the combination 

n(g2) = n(g2) - n*^"(g2) _ n'°^(g2) (11) 

is polynomial in z and 1/z, i.e. in the small and high energy region. Furthermore 
no logarithmic singularities may be introduced at threshold. 



4. Perform a conformal mapping. The change of variables [|T4 



1 + ^ 



(12) 



maps the z plane into the interior of the unit circle of the w plane. Thereby the cut 
[1, 00) is mapped to the perimeter. 
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5. Define [§ 



(13) 



where for n(g^) the terms up to order must be known. The available informa- 
tion transforms into P„(— 1) and Pj^'^O), (A; = 0, 1, . . . , n + no — 1), where uq is the 
number of moments for z 0. 

6. In the last step a Pade approximation is performed for the function Pn{uj). This 
means that Pn{i^) is represented through a function 

r./.n/ X ao + aiu;H h ajtu* 

1 + biu + ■ ■ ■ + bjuj^ 

where the number of coefficients on the r.h.s. depend on the amount of information 
available for Pn{uj). 

Once an approximation for P{uj) is known the above steps have to be inverted in order 
to obtain the function n(g^). 

Due to the structure of Eq. (|1^) some Pade approximants develop poles inside the unit 
circle (|a;| < 1). In general we will discard such results in the following. In some cases, 
however, the pole coincides with a zero of the numerator up to several digits accuracy. 
These Pade approximations will be taken into account in constructing our results. To be 
precise: in addition to the Pade results without any poles inside the unit circle, we will 
use the ones where the poles are accompanied by zeros within a circle of radius 0.01, and 
the distance between the pole and the physically relevant point g^/M^ = 1 is larger than 
0.1. 



3 Current correlator in the effective theory 

The computation of the current correlators is performed in the full theory. In this Section 
the connection to the correlator in the HQET is described. This connection is exploited 
to obtain the leading logarithmic terms of R'"{s) and R^{s) at threshold. 

In the effective theory we denote the MS renormalized current which couples to a 
massive and a massless quark by 

JT = qTQ, (15) 

with r G {1, 7^^, ^75, 7^75}. The light quark flavour q is identical in the full and the 
effective theory which is not true for the heavy one as indicated by the tilde. In the 
effective theory Q is considered as a static quark which fulfills the relation 

Q = 7oQ. (16) 
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Due to this condition the set of 7 matrices can be divided into two groups: one which 
commutes ({^75,7j}, j = 1,2,3) and one which anti-commutes ({1,7^-75}) with 79. 

Let us next consider the corresponding correlators which are defined in analogy to 
Eqs. d) and (|) 



nri,r2(go) 



z /dxe^'^^(0|Tjr,(a;)j^ rO)|0) 



(17) 



They only depend on the zeroth component go as due to the Feynman rules of HQET in 
coordinate space one has {0\Tjr{x)jl{0)\0) ~ 6{x). Furthermore, in our case where only 
one of the quarks is massive, one can write all correlators in terms of n(go) = ni^i(go)- It 
is possible to show that the following equations holdQ 



n 



175,17. 



5(^0) 



n^,.'v,(go) 



n. 



7»75,7j7; 



5(^0) 



n(go) , 



Thus in HQET there is only one independent correlator, n(go), which will be considered 
in the following. The corresponding current will be denoted by j. 

The relation between j and {S = v, s) has been computed in |]T5[ up to order a^. 
For /X = M it can be written in the form 



/ = Cs{M)j, 

where the decouphng constants Cs{M) are given by 



(19) 



a(M) 



("/) 
as 



(M) 



+ 



("/) 
as 



(M) 



2 ^ 



TT 



TT 



/1453 173 



7 



11 



+ CfCa 



6821 21 7, , 

\ C2 C2ln2 H Ca 

2304 16^ 4^ 16^ 



V 768 4. 

CpTrii 



C2 + -C2ln2-— Cs 
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2 

/445 1 
l576 + 4^^ 



+ CfT 



/709 
V576 



(20) 



Cs(M) 



1 + 



+ 



+ CfCa 
+ CfT 



TT 

/1351 



TT 



I- 



768 
/169 



3 3 
— C2 + -C2 In 2 
16^ 4^ 



1 

' 16~ 



369 
256 



15. 3, , 1 \ 
+ — C2 - -C2 n2 - -C3 
16^ 2^^ 8^V 



C3 + CfTui 



95 



192 4 



Vl92 



(21) 



with ^2 = 7r^/6 and (3 ~ 1.202 057. The superscript attached to as defines the number of 
active fiavours. j is still defined with nj active fiavours. The transition to a theory with 

^ This is true only if no power suppressed non-perturbative corrections are taken into account, what 
is imphcitly understood throughout the present paper. 
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only Til active flavours is achieved through 



J = c{M)y 



(22) 



with [116 



89 /^("'HmV 



576 



TT 



(23) 



where again = M has been adopted. In analogy to Eqs. (plTf ) and (|T^) one can define 
polarization functions also in the primed theory. Again there is only one which is inde- 
pendent. It will be denoted by n'(5'o) which only depends on the massless quark flavours. 
Thus besides the renormalization scale there is only one more dimensionful quantity 
which in the framework of HQET is usually chosen to be = ^/s — M. As a consequence 
the renormalization group equation for j' takes the simple form 



d/i' 



-J 



with 17 



7' 



TC 



TT 



-CM — - -C2 



CaCf 



CpTrii — 
48, 



49 



192 8 



(25) 



Note that the primed quantities only depend on a^"'-*. This allows for a simple recon- 
struction of the logarithms ln(yU^/ci;^) at O(a^) of R' deflned through 



R'iuj] 



2nlm 



n'(go 



(26) 



Once they are at hand Eqs. (^), ( pTD and (|^) can be used to predict the logarithms of 
and R^ at threshold via the relations 



R^{s) 
RHs) 



C,{M)C{M)]'^R' + 0{v') 



m{M) 
M 



C,{M)C{M) 



^'^R' + Oiv' 



(27) 



where it is understood that uj is expressed in terms of v via the relation 



1-^ 
M 



(28) 



an expansion for f — > is performed and only the leading term is kept. 
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At this point we want to mention that the developed formahsm — in particular the 
decoupling relations ([19D and (p3D — only applies to the imaginary part. For the polar- 
ization function one would have to take into account additional contributions. 

The procedure described above determines the logarithmic parts of the leading thresh- 
old term. They are incorporated into the Fade procedure as described in Section (cf. 
point 0). Note, that the non-logarithmic part for are not fixed via this procedure. In 
Section ^ we will extract numerical approximations from our Fade results. 

At the end of this Section we want to list the available information for R'{uj). Using 
Eqs. dH) and (H) and the 0{as) result for R' @, || one obtains 



R'iij) 



+ 



TT 
2 r 



TT 



/141 11 



19 



+ CaCf + — - it + — C2 + —L 



V 32 



12 



11 



32 



13 1 



+ CpTrii ( + ( -— + - ln2 - -C2 ) - -L, 



(29) 



with = ln(/i^/ct;^). The coefficients cff, cfa and cfl are not known. In Section |^ we 
will provide numerical approximations. 

With the help of Eq. (p7|) one obtains the leading terms of R^ and R^ for f — > 0. 
Separated according to the colour factors they read {n = M). 



R 



'V,thr 



R 



•s,thr 



TC 



27 



+ I2C2 - 91n2 - 91nw 



+ 



TT 



Cp 1 c^pp + 



147 



27 



27 



3OC2 + — In 2 Int; + — In^ w 



+ C^Cp (4^ + (-^ - 19C2 + y ln2) Int; + ^ In^ 
+ CpTrii [cIl + + 8C2 - 6 ln2) Int; - 3 In^ v 



+ CfT f^-lOC: 



V 



(30) 



TT 



TT 



+ 



+ CaCf e^pA 



C^(13 + 8C2-61n2-61nt;) 



73 



Cf[c'pp+ -— -20C2 + 91n2 lnw + -ln"t; 



141 38 



4 



38 \ 11 \ 

yC2 + llln2j lnt; + y In^ v\ 



8 




Figure 1: Sample diagrams contributing to the current correlator of Eqs. (|T]) and (^. 
The straight and loopy lines represent quarks and gluons, respectively. One of the quarks 
connected to the external current carries mass M whereas the other one is massless. 



+ CFTni[c'pL + (l3 + yC2-41n2)lni;-21n2i; 



+ CfT 



/727 



(31) 



with ^2 = The smooth behaviour proportional to is a consequence from the fact 

that jR! is, for dimensional reasons, proportional to cj^. This is in contrast to the diagonal 
correlators where at order aj. even 1/v singularities appear |]T|. Note that in Eqs. 



and (|3TD the terms proportional to the colour factor CpT are completely fixed. For later 
use it is convenient to display explicitly the relations between c and c^: 



-FF 



--FA 



^FL 



-FF 



-FA 



'-FL 



1427 

Qcff - + 181n2 



269, , „ 33, 

— C2 + 42C2ln2-yC3 



6c 



FA 



64 

6821 63, ^ , ^ 27, 

\ C2 - 2IC2 n2 H Cs 

192 4 ^ 4 ^'^ 



-92.3884 + Qcff ■ 
-25.4483 + 6cfa , 



445 

6cfl + ^ + 2*^2 



+ 14.2056 + 6cfl 



(32) 



4£ff - ^ + 61n2 - + 12C2ln2 - 7(3 

45FA-^ + ^C2-6C2ln2 + ^C3 ~ -8.79653 + 45^^ 
96 2 2 



-24.1011 + 4gi.i. 



47 

45FL + - + 2C2 



5.2482 + AdFL ■ 



(33) 



4 Spectral function in full QCD 

In this Section we discuss the computation of the polarization function in full QCD. 
Explicit results are given for the imaginary parts which constitute physical quantities. In 
a first step the input quantities needed for the Fade procedure are provided. 
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In Fig. ^ typical diagrams are pictured. Altogether roughly 30 contributions have to 
be considered at three-loop order. Although their number is relatively small we used 
GEFICOM [jl8[ for the automatic computation. GEFICOM uses QGRAF Q for the generation 
of the diagrams. In case an asymptotic expansion has to be applied LMP ||20| or EXP |^T[ 



are 



used for the generation of the sub-diagrams. The occurring vacuum diagrams are passed 
to MATAD and the massless propagator type diagrams are evaluated to MINCER . 



More details on the automatic computation of Feynman diagrams can be found in . 

At three-loop order terms up to order (2;^) could be evaluated for the vector 
(scalar) correlator in the limit 2; — 0. For convenience we list below the analytical 
results for the one- and two-loop correlator up to order and present the results for 
the n'-^''''^ in the Appendix. The analytical results can also be found under the URL 
http://www-ttp.physik.uni-karlsruhe.de/Progdata/ttp01/ttp01-14. In the low- 
energy limit we obtain at order a\ and a\ 



n(i).^ 



167r2 
3 

IGtt^ 



1 2 2 1 3 1 4 1 5 2 
-zA z^ A z^H 

2 15 18 35 60 189 



2 6 

z^ + 



1 7 

z 



1 , 25\ 
2^^ + 48j^ + 
1199 \ 4 



12600 
222433 \ 



^ (l40^^^ 7056000 



167r2 
3 

167r2 



1 



15 



+ ... , 

1 



1! 

37883 
604800 



140 
503 



+ ... , 



2 . 5 \ 2 / 1 . 

— C2 H -2 + — C2 + 

^ 18/ V18^ 3240 



z 



z' + 



2 ^302^N 
529200/ 



189 



(34) 



+ -z + -z'^ + — z^ ^ — z" ^ 2° + 



15 



1 6 1 7 

— z^^ z^ 

84 126 



i—C — 
^ V30 ^ 216 



30 ~ ' 105 

3^^+12) ^+(6^^ + 48) ^ +(15^^+72)^ 

1 

8"4 



+ 



12/ V6^ 48 

/ 2 ^ 7001 \ 



' ' 75600) + 



Vl26 423360 



18961 \ 7 



V105 
+ ... , 



38113 
604800 



(35) 



where the on-shell quark mass M has been used as a parameter. 

In the high energy region eight expansion terms could be obtained both in the vector 
and scalar case. Again we provide in the main text the one- and two-loop results in 
analytic form and refer for the analytic three-loop terms to the Appendix and to the URL 
http : / /www-ttp . physik . uni-karlsruhe . de/Progdata/ ttpOl/ttpOl-14. Using again 
the on-shell mass the one-and two-loop results read 



167r2 
1 



16 8 , 

\- -L^ 

9 3 



+ 



10 z^ 



1 1/4 
4L^ --2-+ -L, 
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^Lt-- -+ -L,--L 
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25 45 ^ 81000/ 
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36000. 
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1 2 1 
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16 
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1 



8000/ 



10243 \ 1 
108000/ ^ 
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11025 



4 ^2 532267 \ 1 
105 ^ 9261000/ 



(36) 



(37) 



with L, = -(ln(-z))/2. 

In order to incorporate the available information at threshold one has to perform an 
analytical continuation of the expressions in Eqs. (|30|) and (|3l|). Taking the logarith- 
mic parts of Eqs. (|30|) and (31) one obtains the quadratic and cubic logarithms for the 



polarization functions which read 
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'-'-FH,log 



3 3 



167r2 2 
3 

167r2 
3 

167r2 
3 

167r2 
0, 



1 



^ ( 1- - 

z 
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Note that, except for the coefficient of the quadratic logarithm in the structure FF, there 
is agreement between the coefficients of the vector and scalar correlator. The subscript 
"/og'" reminds that the imaginary parts of the expressions in Eq. (|38|) only reproduce the 
logarithmic parts of ( PUj ) and (^Tj). We want to mention that the polarization function 
itself also contains linear logarithmic terms which are not yet known at order a^. Below 
we will derive numerical estimates for them. For the two-loop correlators the linear 
logarithms are known. They can be obtained from the constant contributions in Eqs. (|30|) 
and (|T|). 

Now the complete input is available and the individual steps described in Section ^ 
can be performed leading to a large variety of Fade approximants. For the results we 
present in the following those Fade approximants are chosen which contain for their 
construction at least terms of order z^ and 1/z^ in the small and large momentum region, 
respectively. Furthermore we demand that the difference of the degree in the numerator 
and denominator in Eq. (|^) is less or equal to two. 

In Figs. and ^ the two-loop results are plotted as a function of v and x, respectively. 
In addition to the roughly 15 Fade results also the exact expression is plotted. However, it 
is not possible to detect any difference — even close to the threshold at f = 0. The dashed 
lines in Fig. |^ correspond to the results of the high-energy expansion including terms up 
to order 1/^''. Although in these curves only the information from x ^ is incorporated 



12 




Figure 3: R^^^''"{s) and R^^^''^{s) as a function of x. 




Figure 4: Rpp{s), Rp/'{s), Rpi{s) and Rpjf{s) as a function of v. 



one observes a perfect agreement with the exact results up to x ~ 0.9 {v ~ 0.10). 

The Fade results for the individual colour structures at order are plotted in 
Figs. 1^, P and 0. Again, despite the fact that in each plot approximately 15 curves are 
shown no difference between them can be observed. In those plots where x is chosen for 
the abscissa also the result of the high-energy expansion containing terms up to order 
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Figure as a function of x. 

is plotted as a dashed line. There is excellent agreement with the semi-numerical result 
at least up to x ~ 0.5, which corresponds to f ~ 0.60. In some cases even an agreement 
up to a; ~ 0.7 is observed. This corresponds to f ~ 0.34 which is already quite close to 
the threshold. 

A remark concerning the plots for Rp\'^ and R^^^'^ are in order. In Fig. ^ it can be 
seen that for a; — they tend to +oo and —oo, respectively. However, the turn-over takes 
place in a small region of x which is beyond the resolution in Fig. ^ where the variable v 
is used and thus the high-energy region, i.e. the region for t> — ^ 1, gets squeezed. 

The analytical formulae which result from the semi-numerical Fade procedure are 
quite long. Thus we refrain from listing them explicitly. Instead, a typical representative 
for the two-loop results and for each colour structure at three-loops can be found under 
the URL http : / /www-ttp . physik . uni-karlsruhe . de/Progdata/ ttpOO/ ttpOO-25. 

5 Spectral function in HQET 

In this Section we concentrate on the behaviour close to the threshold. In particular we 
want to extract the leading non-logarithmic term which is of order f ^ and afterwards 
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transform the result to HQET. 
The expressions in Eq. (|3| 



are constructed in such a way that the combination 
Im[n(g^) — Ilj^g{q'^)]/v'^ approaches a constant for v —>■ 0. However, in the method 
we use for the computation of n(g^) it is not possible to incorporate the constraint 
that Im[n(g2)] ~ v"^ in analytical form. Our procedure provides rather a function 
which numerically imitates the quadratic dependence. As a consequence the expression 
Im[n(g^) — nf^^(g^)]/f ^ diverges for very small velocities and provides quite different nu- 
merical values for the different Fade approximants. Nevertheless one could try to extract 
numerical approximations for the non-logarithmic coefficient. 
For this aim we consider the quantity 



nslm[Ilris + ^e)]-R^iff;^i 



s=M'2(l+v)/{l-v) 



where = Un, n, = Svr and X e {F, FF, FA, FL, FH}. R^^ff^"^ is obtained from 

Eqs. (|30D and ( ^I]) by setting the unknown constants to zero. We want to use t'^'^{v) 
in order to fit the coefficient in Eq. (|39D for each individual Fade approximant. This 



(39) 
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Figure 7: Rpp^{s), Rp]l'^{s), Rpp^{s) and Rp]f{s) as a function of x. 

requires small values of v. However, due the very structure of the Fade result the ratio 
Tx^'^{v)/v'^ becomes unstable for f — 0. Thus, in a first step we determine a minimal value 
Vq for which T^'^ {v) /v"^ evaluated for all Fade results still show reasonable agreement. 
Furthermore we determine an effective coefficient from the equation 

T«''(t;o) = iV,4,,ff^;o^ (40) 

and require that the error on c^xes froin the different Fade approximants is less than 10%. 
The results for vq and c^eff f*^^ each colour factor both for the vector and scalar current 
can be found in Tab. |I[ 

In order to account also for higher order terms we choose in a next step a value Vi > vq 
in such a way that, first, f i — f o ~ fo and, second, Tx^'^ {v) /v^ displays an approximately 
linear behaviour for Vq < v < Vi. Then a two-parametric fit is performed in the interval 
[fct'i], i.e. the coefficients and d^x determined from the equation 

tF(^) = Kv^ [c'x + d'xv) . (41) 

Finally Eqs. (p2D and (^3]) are used to convert the result for to a numerical estimate for 
cx- The corresponding results for the individual structures can again be found in Tab. [|. 
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Table 1: Results from the fits for the coefficients as described in the text separated accord- 
ing to the colour factors and the tensorial structure of the current correlator. The errors 
indicated in the brackets arise from the comparison of the different Fade approximants. 
They are omitted in the case of cx as they are much smaller than the systematic error of 
the extraction procedure. 

The corresponding errors are omitted in the case of cx since they are much smaller than 
the systematic error of the extraction procedure. 

In order to get confidence in our prescription for the determination of let us have a 
look at the two-loop results. The values given in Tab. |^ can be compared with the exact 
results which read 

cl = 27.00, 
4 = 22.00, 

cf = 6.50. (42) 

One observes a very good agreement of the approximated results for both cp and Cp 
with the exact ones0. Although the error induced by the comparison of the different 
Fade approximations is slightly smaller than the deviation from the exact result the 
considerations at order as are quite promising for the determination of the coefficients at 
order af. 

As was already mentioned above the coefficients completely determined. Thus, 

also they can be used in order to test our prescription. However, both in the vector and 
scalar close to zero which, of course, leads to large relative errors. Nevertheless, 

let us present the results which read for vq = 0.05 and vi = 0.10 

= 0.20(6), 
4^ = 0.16(10), 
4^,eff = 0.50(6), 

c'fh = 0.44(11). 

^This is also the case for vq = 0.03 and vi — 0.06. 

•^As can be seen in Figs. I and I both R'^plf and i?^^]^" are rather smooth around v 
for larger values of vq and vi . 
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Table 2: Results from the fits for the coefficients as described in the text separated for 
different values of rii. As far as the errors are concerned the same statements hold as in 
Tab. I. 



The comparison with the exact values obtained from Eqs. (^) and ( |3lD 

c^H = 0.175... , 

4^ = 0.455 . . . , (44) 

shows that the agreement of the central values is within 10%. 

In order to obtain our final predictions for the coefficients cx we proceed as follows: 
the coefficients cx given in Tab. ^ have to be independent from the tensorial structure of 
the current correlator from which they are determined. On the other hand, guided by the 
numbers shown in Tab. 0, we can define criterions concerning the stability and reliability 
of the extraction procedure. As a first criterion we require that Vq should be as small as 
possible. This clearly suppresses higher order terms in v. Furthermore we consider the 
ratio dx/cx and regard that colour structure as more reliable where the ratio is smallest. 
As dx represents an effective constant containing the effects of the higher orders in v also 
this criterion selects the structure where they are suppressed as compared to the f ^ terms. 
Following these rules we choose the vector correlator in order to obtain cpA and the scalar 
one for cpp and cpL- This leads to the following results 

cfp = 21(6), 
Cfa = 1.2(4), 

CFL = -2.3(7), (45) 

where we assigned a conservative error of 30%. Note that all numbers given in the last 
column of Tab. |I| are consistent with our final predictions of Eqs. (^S|). 

In practical applications it is often sufficient to know the result at order for the 
numerical value A^^^ = 3 which motivates the following procedure: in a first step the 
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moments of the colour structures FF, FA and FL are added for fixed rii and the Fade 
procedure is performed. Afterwards the prescription described around Eq. (^) is apphed. 
The results of the corresponding analysis can be found in Tab. ^ which allow to deduce 
the following compact expression for c„j which governs the dependence on nf. 



-■ni 



46(15) - 1.2(4) n/, 



(46) 



At first sight the errors appear quite large. However, one has to recall that only the large- 
and small-g^ behaviour of the polarization function serves as input whereas the quantity 
in Eq. (|6) corresponds to the imaginary part at = M^. Note also that at order 
the numerical values for vq have to be chosen larger than at order which increases the 
influence of the higher order terms in v. 

To summarize: the complicated threshold structure of the Fade approximants allows 
only a rather rough determination of the constants and c^. However, the peculiar 
structure of Eqs. (|3^) and ( p3D connecting to c makes it possible to extract the universal 
coefficient c with a reasonable accuracy. 

Finally we are in the position to write down the expression for B!{uj) up to order a^. 
Inserting the colour factors into Eq. (|29|) and using the results from Eq. (^61) one obtains 



1 + (8.667 + L^) + [ -^-^ (46(15) + 35.54L, 



7r 



TT 



+ 1.875L2 + m (-1.2(4) - 1.583L^ - 0.08333L^ 



(47) 



We can conclude that the coefficient at order ai is quite large and has a mild dependence 
on ni. 



6 Applications and conclusions 

As an application of the vector and axial-vector current correlator we want to discuss the 
single-top-quark production via the process qq tb. Some sample diagrams contributing 
to this process are shown in Fig. ^. 

The corrections of order as to the (total) single-top-quark production rate are quite 
large. They amount to about 54% and 50% for Tevatron and LHC energies, respec- 
tively ||25[, where 18%, respectively, 17% arise from the final state corrections. This 
makes it necessary to consider also the corrections of order a^. Due to the appearance 
of an interference between the initial and final state (c.f. Fig. §(f)), which for the first 
time happens at order a^, the complete calculation is a non-trivial task. However, with 
the results of this paper we are in the position to perform a first step and consider the 
leading term in the large- iVc expansion. 

One observes that the contributions of the diagrams where gluons connect the initial 
and final states are suppressed by at least a factor l/N^ in the large Nc limit as compared 
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Figure 8: Sample diagrams contributing to the process qq 
lines represent W bosons and gluons, respectively. 



th. The wavy and loopy 



to the diagrams in Fig. ^(c)-(e). For the latter, together with the contributions of Fig. 
and (b), the differential cross section can be written in factorized form 



da 

dg' 



2 {PP 



tb + X) = a{pp^W* +X) 



imn 



w 



{q\MlMl 



Ti{q^ 



(45 



where Uw corresponds to the transversal part of the W boson self energy which is con- 
nected to the vector correlator of Eq. (|I|) through 

Uw{q') = V2GFM',,\V,,\yW{q'). (49) 

At order there are also diagrams like the one in Fig. H(g) which appear for the 
first time. In principle they also lead to the same final state as the diagram in Fig. §(e). 
However, one has to note that the W boson generating the top and bottom quark is 
radiated from a light quark flavour. This suggests that their contribution is small although 
there is only a suppression by a factor l/N^ as compared to Fig. ||(c) and not by 
like for the diagram in Fig. Kf). 

Thus, if we restrict ourselves to the leading term in it is possible to use the results 
for K" obtained above in combination with Eq. (|^) to perform a theoretical analysis at 
order to the single-top-quark production in the large- A^c limit. In order to obtain the 
total cross section the corresponding parton distribution functions would be needed to 
the same order. 

The production cross section of the virtual W* boson is identical to that of the Drell- 
Yan process qq — > eug- The latter is known to O^a^) from Ref. Thus we can take the 
proper ratios to make predictions in the large- iV^ limit at NNLO free from any dependence 
on parton distribution functions. As an example, we can consider the following expression 

{pp ^W* ^ th) 



da 



_ Im [n,,(g^)] 

W* -> ez/e) Im [ne.(g2)] 

= NAVtt\''R''(s) 



(50) 



In Fig. 0the LO, NLO and NNLO results ofR'is) are plotted in the range = 200 . . . 400 
GeV where Mt = 175 GeV and as{Mz) = 0.118 has been chosen for the numerical anal- 
ysis. Whereas the 0{as) corrections are significant there is only a moderate contribution 
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Figure 9: LO (dotted), NLO (dashed) and NNLO (solid) results of i?^(s). 



from the order terms. In the range in shown in Fig. ^ they are below 1% of the 
Born result. Note that the NNLO correction to the Drell-Yan process are also small and 
amount to at most a few percent (see e.g. Thus, in case there is no kinematical 

magnification for the diagrams in Fig. ^(f) and (g) we can conclude that the radiative 
corrections to the single-top-quark production via the process qq — *■ fb are well under 
control. 

The scalar and pseudo-scalar correlator covers properties connected to a charged Higgs 
boson. The latter appear in theories beyond the SM which are usually characterized by 
an enlarged Higgs sector containing Higgs bosons with different quantum numbers. For 
example, one of the most appealing extensions of the SM, the Minimal Supersymmetric 
Standard Model (MSSM), contains two complex iso-doublets with opposite hyper-charge 



(see e.g. p8[), resulting in five mass eigenstates of (pseudo-) scalar physical Higgs fields: 
two neutral CP-even [H^ and /i"), one neutral CP-odd {A) and two charged {H^) Higgs 
bosons. 

Let us consider a generic charged Higgs boson coupled to fermions through 



H+DU 



{V2G1 



1/2 



(51) 
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where the corresponding quark current is given by 

Jm = ^^[a(l-75) + 6(l + 75)]^. (52) 
Here U and D represent generic up- and down-type quarks, respectively, with MS masses 



mu and mo = 0. Eq. (|5^) only covers the contributions from a boson; the formulae 
for a Higgs boson with negative charge are analogous. The parameters a and h are model 
dependent and are left unspecified. 

The decay rate of the boson into quarks and gluons can be written in the form 

V{H+^UD) = v^CpMH+Im [njy(M^+)] , (53) 

where Mu is the pole quark mass and 11// (g^) is given by 

q'UHiq') = Jdxe''^^TJH^{x)JH-m = {a' + b')q'U'{q'), (54) 

Thus, we arrive at the following expression for the hadronic decay rate of the charged 
Higgs boson 

T{H+^UD) = ^^^MH+{a^ + b^)R'{Mfj^). (55) 



In Fig. |T0] W{Mjj+) is plotted at LO, NLO and NNLO. Again it turns out that the 
radiative corrections are well under control as order terms contribute at most of the 
order of 1%. 

As an application of the correlator in the effective theory we want to mention the 
determination of the meson decay constants via QCD sum rules where the Borel transform 
of B! as given in Eq. ( ^TD enters as a building block. Besides the perturbative part the 
sum rules also obtain contributions from non-perturbative condensates which, however, 
are numerical less important The typical scale which has to be used in Eq. ( ^Tj ) is of 
the order of 1 GeV |^ which leads to sizeable corrections both at order and at order 

To be more precise let us choose af\Mz) = 0.118 which, using two- loop accuracy, 
leads to q;^^^(1.3 GeV) ~ 0.37 For u = 1.3 GeV and n/ = 4 the order as corrections in 
Eq. ( |i7| ) amount to about 100%. The terms at 0{a'^) contribute with additional 60(20)% 
where the sign is the same as for the LO correction. A careful analysis is necessary in 
order to decide whether perturbation theory can safely be applied in this case. 

To conclude, in this paper the non-diagonal current correlator formed by a massive and 
massless quark has been considered. Moments in the low and high energy region have been 
computed analytically in full QCD. Furthermore the leading logarithmic contributions 
arising at the quark threshold q^ = have been obtained from the reconstruction of 
the logarithmic terms of the spectral function in the effective theory. This information is 
combined with the help of conformal mapping and Fade approximation to obtain semi- 
numerical results for the vector and scalar correlator and in particular their imaginary 
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Figure 10: LO (dotted), NLO (dashed) and NNLO (solid) results of R'{s), M = Mt = 
175 GeV. 

parts valid for all values of q and M. In a next step the various Fade results are used 
in order to obtain the leading non-logarithmic coefficient at threshold which is in turn 
transformed to the effective theory leading to a prediction of the spectral function up to 
order (see Eq. (|47|)). As applications we considered the effect of the correction to 
the single-top quark production and the decay of a charged Higgs boson. 
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Appendix: 

Analytical results for the moments at order a 



In the limit 2; — > the results for n^^^''' and n^^^'* parameterized in terms of the on-shell 
mass read 
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where = JVP has been adopted. In the high-energy hmit the corresponding results read 
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In the above expressions the following symbols have been used: ^4 = Li4(l/2), 
Ci = Cl2(7r/3) = Im[Li2(exp(m/3))], C2 = ttVG, C3 ~ 1.202 057, C4 = ttVQO, (5 ~ 
1.036 928 and = — (ln(— z))/2. A Mathematica input can be found under the URL 
http : //www-ttp . physik . uni-karlsruhe . de/Progdata/ttpOl/ttpOl-14. 
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